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Abstract. The rehc abundances of the hght elements synthesized during the first 
few minutes of the evolution of the Universe provide unique probes of cosmology and 
the building blocks for stellar and galactic chemical evolution, while also enabling 
constraints on the baryon (nucleon) density and on models of particle physics beyond 
the standard model. Recent WMAP analyses of the CBR temperature fluctuation 
spectrum, combined with other, relevant, observational data, has yielded very tight 
constraints on the baryon density, permitting a detailed, quantitative confrontation of 
the predictions of Big Bang Nucleosynthesis with the post-BBN abundances inferred 
from observational data. The current status of this comparison is presented, with an 
emphasis on the challenges to astronomy, astrophysics, particle physics, and cosmology 
it identifies. 



1 Introduction and Overview 

Our Universe is observed to be expanding and filled with radiation (the Cosmic 
Background Radiation: CBR), along with "ordinary" matter (baryons = nucle- 
ons). It is well known that if this evolution is followed backwards in time, then 
there was an epoch during its early evolution when the Universe was a "Primor- 
dial Nuclear Reactor", synthesizing in astrophysically interesting abundances 
the light nuclides D, '^He, '^He, and ^Li. Discussion of BBN can start when the 
Universe was some tens of milliseconds old and the temperature (thermal en- 
ergies) was of order a few MeV. At that time there were no complex nuclei, 
only neutrons and protons. Since there are nearly ten orders of magnitude more 
CBR photons in the Universe than nucleons, photodissociation ensures that at 
such high temperatures the abundances of complex nuclei are negligibly small. 
However, as the Universe expands (and the weak interactions transmute neu- 
trons and protons), collisions among nucleons begin to build the light nuclides 
when the temperature drops below ~ 100 MeV, and the Universe is a couple of 
minutes old. Very quickly, almost all neutrons available are incorporated in the 
most tightly bound of the light nuclides, '*He. As a result, the ^He primordial 
abundance (mass fraction: Yp) is a sensitive probe of the competition between 
the weak interaction rates and the universal expansion rate (the Hubble param- 
eter: H); ''He is a cosmological chronometer. The reactions building ^He are 
not rate (nuclear reaction rate) limited and, therefore, Yp is only weakly (log- 
arithmically) sensitive to the baryon density. In contrast, the BBN abundances 
of the other light nuclides (D, ^He, ^Li) are rate limited and these do depend 
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sensitively (to lesser or greater degrees) on the baryon density; D, '^He, and ^Li 
are all potential baryometers. 

The relic abundances of the light nuclides predicted by BBN in the "stan- 
dard" model of cosmology (SBBN) depend only on one free parameter: the nu- 
cleon density. There are, therefore, two complementary approaches to testing 
SBBN. On the one hand, the primordial abundances inferred from observational 
data should be consistent with the SBBN predictions for a unique value (or 
range) of the nucleon density. On the other hand, if there is a non-BBN con- 
straint on the range allowed for the baryon density, this should lead to SBBN- 
predicted abundances in agreement with the observational data. Is this the case? 
If not, why not? That is, if there are conflicts between predictions and obser- 
vations, is the "blame" to be laid at the feet of the observers (inaccurate data 
and/or unidentified systematic errors?), or the astrophysicists (incorrect mod- 
els for analysing the data and/or extrapolating from abundances determined at 
present ("here and now") to their primordial ("there and then") values), or are 
the standard models of particle physics and/or cosmology in need of revision? 

1.1 The Status Quo: Observations Confront SBBN 

Space limitations prevent my presenting a full-fledged review of the history of 
the observational programs along with the evolution of the comparisons between 
theory and data. For some recent reviews of mine, see ^ and further references 
therein. Instead, an overview is presented which highlights the challenges to 
SBBN. The remainder of this article is devoted to some of the key issues associ- 
ated with each of the light nuclide relic abundances. 



1 


1 


WMAF 




'Li 














D 


•He 

1 






2 













Fig. 1. The baryon density parameter, 7710, inferred from SBBN and the relic abun- 
dances of D, ^He, ^He, and ^Li (filled circles), along with the non-BBN determination 
from WMAP (filled triangle). See the text for details. 
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If the standard model (SBBN) is the correct choice and if the primordial 
abundances inferred from the data were free of systematic errors, then the baryon 
densities determined from D, "^He, ^He, and ''Li should agree among themselves 
and with that inferred from the CBR (WMAP) and other, non-BBN, cosmo- 
logical data. In Figure 1 are plotted the various values of the baryon density 
parameter^ determined by SBBN and the adopted primordial abundances and, 
also, from the WMAP-team analysis As may be seen from Fig. 1, the SBBN 
D abundance is in excellent agreement with the WMAP-inferred baryon den- 
sity. However, neither ^Li nor ^He agree with them or, even with each other. 
While ^He is consistent with SBBN deuterium and with WMAP (and is not in 
disagreement with ''Li), the very large uncertainty in its primordial abundance, 
combined with its relative insensitivity to the baryon density, render it - at 
present - an even less sensitive baryometer than is ''He. 

In the next section each light nuclide is considered in turn, its post-BBN 
evolution briefly reviewed along with identification of a few of the potential 
challenges to accurately inferring the primordial abundances from the observa- 
tional data. Then, having established that the current data - taken at face value 
- are not entirely consistent with SBBN, I investigate whether changes in the 
early universe expansion rate can reconcile them. 

2 Primordial Abundances — Evolution, Uncertainties, 
Systematics 

While Figure 1 suggests some problems with SBBN, the optimist may prefer to 
conclude that observations have provided impressive confirmation of the stan- 
dard cosmological model. After all, if the standard model - or something very 
much like it - were not correct, there'd be no good reason why the abundances 
of four light nuclides, which range over some nine orders of magnitude, should 
be just such that the baryon densities inferred from each of them lie within a 
factor of three of each other, in nearly perfect agreement with that derived inde- 
pendently from non-BBN data. Only recently, when cosmology has entered an 
era of great precision, has it become important to distinguish accuracy from pre- 
cision and to revisit the path from precise astronomical observations to accurate 
abundances. For each of the light nuclides of interest here, an all too abbreviated 
overview of the current uncertainties is presented, with the intentional goal of 
creating controversey in order to stimulate future observations and theoretical 
analyses. 

2.1 Deuterium 

Deuterium is the baryometer of choice. During its post-BBN evolution, as gas is 
cycled through stars, D is only destroyed (setting aside rare astronomical events 

^ After annihilation during the early evolution of the Universe, the ratio of baryons 
to photons is, to a very good approximation, preserved down to the present. The 
baryon density parameter is defined to be this ratio (at present): = n^/n^; ryio = 
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where, far from equilibrium, tiny amounts of D may be synthesized). Therefore, 
observations of D anywhere, at any time {e.g., the solar system or the local 
ISM), provide a lower bound to its primordial abundance. As a result, it is ex- 
pected that observations of D in regions which have experienced minimal stellar 
evolution {e.g., the high redshift and/or low metallicity QSO Absorption Line 
Systems: QSOALS) should provide a good estimate of the relic abundance of deu- 
terium. Kirkman et al. |3] have gathered the extant data; see |3] for details and 
related references. In Figure 2 are shown the QSOALS deuterium abundances 
as a function of metallicity; for reference, solar system and ISM D abundances 
are also shown. 

While the observers are to be commended for their heroic work in identifying 
and analysing the tiny fraction of QSOALS which can be used to infer a low- 
metallicity, high redshift, D abundance, there are several unsettling aspects of 
the data displayed in Figure 1. Perhaps most noticeable is the paucity of data 
points. Without at all minimizing the difficulty of finding and analysing such 
systems, it is very nearly a sin to claim that the primordial abundance of a 
cosmologically key light nuclide is determined by five data points. If, however, 
the data points were in agreement within the estimated statistical errors, this 
might be less disturbing. It is clear from Fig. 1 that this is not the case. For 
these five data points the about the mean is ^ 16, suggesting either that 
the uncertainties have been underestimated, or that some of these data may 
be contaminated by unidentified systematic errors. While the dispersion may 
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Fig. 2. The deuterium abundance (by number relative to hydrogen), j/d = 10^(D/H), 
derived from high redshift, low metallicity QSOALS (filled circles). The metallicity 
is on a log scale relative to solar; depending on the line-of-sight, X may be oxygen 
or silicon. Also shown is the solar system abundance (filled triangle) and that from 
observations of the local ISM (filled square). 
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simply be due to the small number of data points, it might be significant that 
the three QSOALS with the lower D/H ratios are LLS, while the two highest 
D/H determinations are fi'om DLAs. 

In the absence of evidence for changing or eliminating any of the current 
D abundance determinations, it is not unreasonable to follow the advice of 
and adopt the weighted mean as an estimate of the primordial D abundance. 
Here, too, there is a (minor) problem. Kirkman et al. |3] advocate finding the 
mean of log(D/H) to determine ?/d = 105(D/H); yo = 10(5+<i°g(D/H)» ^ 2.78. 
In contrast, if the weighted mean of the five D/H determinations is used, j/d = 
2.60. While this difference is well within the dispersion, it reflects the fact that in 
determining the mean of log(D/H), PKS1937 with j/d = 3.24 dominates, while 
for the mean of D/H, HS0105 with ?/d — 2.54 dominates. In what follows I adopt 
j/D = 2.6±0.4, where the error, following is derived from the dispersion in 
D/H determinations. The corresponding BBN (SBBN) prediction for the baryon 
abundance, 7710(D) w G.lj^pg, is shown in Figure 1. 

2.2 Helium-3 

In my talk at this meeting I actually avoided discussion of ^He, until it was forced 
upon me during the question session. In part, this was because this subject was 
ably covered in Tom Bania's talk, in Dana Balser's poster, and in Bob Rood's 
conference summary. In part, however, it was because, in my opinion, for both 
observational and theoretical reasons '^He has more to teach us about stellar and 
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Fig. 3. The '^He abundances (by number relative to hydrogen), 1/3 = 10^(^He/H), 
derived from Galactic Hii regions as a function of galactocentric distance (filled 
circles). Also shown for comparison is the solar system abundance (solar symbol). The 
open circles are the oxygen abundances for the same Hii regions (and for the Sun). 



6 Gary Steigman 



Galactic evolution than about BBN. ■^He is a less sensitive baryometer than is 
D since (D/H)bbn t]^^'^, while (^He/H)BBN k Even more important 

is that, in contrast to the expected, monotonic, post-BBN evolution of D, the 
post-BBN evolution of '^He is quite complicated, with competition between stel- 
lar production, destruction, and survival. For years, indeed decades, it had been 
anticipated that net stellar production would "win" and the abundance of ^He 
would increase with time (and with metaUicity); see, e.g., [Sj. Unfortunately, ob- 
servations of "^He are restricted to the solar system and the Galaxy. Nonetheless, 
since there is a clear galactic gradient of metaUicity (see Fig.|31l, a gradient in ^He 
abundance would also be expected. If net production "wins" , then '^He should be 
highest in the inner galaxy; the opposite if net destruction dominates. The ab- 
sence of any statistically significant gradient in the Bania, Rood, Balser (BRB) 
data 0] (see Fig. points to an extremely delicate balance (cancellation) be- 
tween production and destruction. This suggests that the mean '^He abundance 
in the Galaxy (7/3 w 1.9) might provide a reasonable estimate of the primordial 
abundance. However, BRB recommend that the '^He abundance determined in 
the most distant (from the Galactic center) , most metal poor Galactic H 11 region 
yields an upper limit, t/3 ^ 1.1 ±0.2, to the primordial abundance. The estimate 
of Tjioi^He) » 5. 6^^'° shown in Figure 1 is based on this choice. Had I adopted 
the mean value of 2/3 = 1.9, I would have inferred 7710 (^He) « 2.3. While the 
former choice is in excellent agreement with deuterium (and with ''Li and with 
the WMAP result (2), the very large uncertainty renders '^He an insensitive 
baryometer; the latter option would be consistent with ^He, but not with D (or 
with "^Li or WMAP). 

2.3 Helium-4 

Helium-4 is the textbook example of a relic nuclide whose abundance is known 
precisely but, likely, inaccurately. To be of value in testing SBBN as well as in 
constraining non-standard models, Yp should be determined to 0.001 or better. 
The largest uncertainty in the SBBN prediction is from the very small error in 
the neutron lifetime (r„ = 885.7 ±0.7 s). For the WMAP estimate of the baryon 
density, including its uncertainty, the SBBN-predicted primordial abundance 
is Yp = 0.2482 ± 0.0007, as shown in Figure ^ Also shown in Figure El is a 
record of Yp determinations, from observations of extragalactic, low metaUicity, 
H II regions, covering the period from the late 1970s to the present (2004). During 
this time it has been well known, but often ignored, that there are a variety of 
systematic uncertainties which might dominate the Yp determinations. In the 
hope of accounting for these systematic errors (rather than constraining them 
by observations), the error estimates for Yp have often been inflated beyond 
the purely statsitical uncertainties. Thus, until the late 1990s, the typical error 
estimate for Yp was 0.005. For example, summarizing the status as of 2000, 
Olive, Steigman, and Walker suggested that the data available at that time 
were consistent with Yp = 0.238±0.005; see Fig.0] However, if ^He were used as 
an SBBN baryometer (not recommended!), the error in the baryon density would 
have been ~ 50%. Of course, as the number of Hii regions observed increased, 
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largely due to the work of Izotov & Thuan 6 , the statistical errors decreased. 
For example, from observations of 82 extragalactic Hii regions, in their 2004 
paper Izotov & Thuan quote [H] Yp = 0.2429 ± 0.0009; this data point is shown 
in Fig. 0] In contrast, a very recent, detailed study of the effects of some of 
the identified systematic uncertainties by Olive & Skillman suggests the true 
errors are likely larger than this, by at least an order of magnitude. 

As may be seen from Fig. 0] none of the Yp estimates agree with the SBBN 
prediction, all being low by roughly 2-a. Indeed, from their sample of 82 data 
points Izotov & Thuan derive such a small uncertainty that their central value 
is low by nearly 6-a. In their analysis, Izotov & Thuan commit the cardinal sin 
of examining their data and then, a posteriori, choosing a subsample of 7 H ii re- 
gions to derive their favored estimate of Yp = 0.2421±0.0021. One wonders what 
they may have found from a random series of choosing 7 of 82 data points. In 
any case, this estimate also falls short of the SBBN prediction by nearly 3-ct. Us- 
ing this suspect subsample, Olive & Skillman ^ do find consistency with SBBN 
once they have corrected for the systematic errors they've chosen to include. 
However, one they have ignored, the ionization correction factor almost cer- 
tainly would have the effect of reducing their central value and increasing their 
error estimate. For the entire Izotov & Thuan sample, Olive & Skillman find a 
very large range for Yp, from 0.232 to 0.258 (or, Yp « 0.245 ± 0.013, entirely 
consistent with SBBN)). If the average correction for ionization suggested by 
Gruenwald, Steigman, and Viegas ^ for the somewhat smaller 1998 Izotov & 
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Fig. 4. A summary of the time evolution of primordial "^He abundance determinations 
(mass fraction Yp) from observations of metal-poor, extragalactic Hii regions (see the 
text for references). The solid horizontal line is the SBBN- predicted '*He abundance 
expected for the WMAP (and/or D) inferred baryon density. The two dashed lines 
show the la uncertainty in this prediction. 
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Thuan data set is applied to their 2004 compilation, it would suggest the Olive & 
Skillman central value be reduced and their error increased: Yp « 0.239 ±0.015. 

If, indeed, the true uncertainty in Yp is really this large, it opens the pos- 
sibility that there might be alternate approaches to Yp which are competitive 
with the extragalactic, low-metallicity H ii regions and, even more important, 
complementary in that they are subject to completely independent sets of sys- 
tematic errors. One such example, with a venerable history of its own, is to use 
effect of the initial helium abundance on the evolution of low-mass Popii stars, 
employing the R-parameter (10) . Recently, Cassisi, Salaris, and Irwin jfllj have 
attempted this using observations of a large sample of Galactic Globular Clus- 
ters (GGC) and new stellar models. While they claim an extraordinarily accurate 
determination of Yp (0.243 ± 0.006), this does not seem to be supported by the 
data they present since, for the lowest metallicity GGCs, Y ranges from ^ 0.19 
to ^ 0.27. Nonetheless, if there is the possibility that the R-parameter method 
might achieve theoretical and observational uncertainties ^ 0.01, it is certainly 
an approach worth pursuing. 

An alternate approach, subject to large theoretical uncertainties, would be 
to attempt to use chemical evolution models, tied to the solar helium and heavy 
element abundances, to extrapolate back to the relic abundance of "^He. While 
at present this approach appears to be limited by the theoretical uncertainties 
{e.g., metallicity dependent stellar winds and stellar yields), the following ex- 
ample may serve as a stimulus (or challenge) to those who might believe they 
can do better. In a recent paper employing new yields, Chiappini, Matteucci, 
and Meynet [El find Z\Y=Yq-Yp « 0.018 ± 0.006. Using the recent Bahcall 
& Pinsonneault jE] estimate of Yq « 0.270 ± 0.004, this leads to a primor- 
dial estimate of Yp « 0.252 ± 0.007. Although this result is consistent with the 
SBBN prediction, it would be entirely premature to declare victory on the basis 
of such a crude estimate. The possible lesson illustrated by this example is that 
the error associated with such an approach might not be uncompetitive with 
those from the standard Hii region analyses. 

2.4 Lithium-7 

As with "'He, the recent history of the comparison between the SBBN predictions 
and the observational data leading to the relic abundance of ^Li is one of conflict, 
with the spectre of systematic errors looming large. ^Li, along with ^Li, ^Be, ^°B, 
and ^^B, is produced in the Galaxy by cosmic ray spallation/fusion reactions. 
Furthermore, observations of super-lithium rich red giants provide evidence that 
(some) stars are net producers of lithium. Therefore, to probe the BBN yield of 
^Li, it is necessary to restrict attention to the most metal-poor halo stars (the 
"Spite plateau"). Using a specially selected data set of the lowest metallicity 
halo stars, Ryan et al. |14| claim evidence for a 0.3 dex increase in the lithium 
abundance ([Li] = 12-hlog(Li/II)) for -3.5 < [Fe/H] < -1, and they derive 
a primordial abundance of [Li]p w 2.0 — 2.1. This value is low compared to 
the estimate of Thorburn (1994) [^j, who foimd [Li]p w 2.25 ± 0.10. In the 
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steps from the observed equivalent widths to the derived abundances, the stellar 
temperature plays a key role. When using the infrared flux method effective 
temperatures, studies of halo and Galactic Globular Cluster stars JHI suggest a 
higher abundance: [Li]p = 2.24 ± 0.01. Very recently, Melendez & Ramirez [T7) 
have reanalyzed 62 halo dwarfs using an improved infrared flux method effective 
temperature scale. They fail to find the [Li] vs. [Fe/H] gradient claimed by Ryan 
et al. 2^ and they confirm the higher lithium abundance, finding [Li]p = 2.37 ± 
0.05. As shown in Figure 1, if this were the true primordial abundance of ^Li, 
then ?7io(Li) = 4.5 ±0.4. All of these observational determinations of primordial 
lithium are significantly lower than the SBBN expectation of [Li]p = 2.65^Q'Qg 
for the WMAP baryon density. 

As with "'He, the culprit may be the astrophysics rather than the cosmology. 
Since the low metallicity, dwarf, halo stars used to constrain primordial lithium 
are the oldest in the Galaxy, they have had the most time to modify (by dilu- 
tion and/or destruction) their surface abundances. While mixing of the stellar 
surface material with the interior would destroy or dilute the prestellar lithium 
abundance, the very small dispersion in [Li] among the low metallicity halo stars 
suggests this effect may not be large enough to bridge the « 0.3 dex gap be- 
tween the observed and WMAP/SBBN-predicted abundances ([Li]f,''* w 2.37 
versus [Li]p'^'^'^ « 2.65); see, e.g., |18j and further references therein. 

3 Non-standard BBN 

The path from acquiring observational data to deriving primordial abundances 
is long and complex and littered with pitfalls. While the predicted and observed 
relic abundances are in rough qualitative agreement, at present there exist some 
potential discrepancies. These may well be due to the data, the data analysis, 
or the extrapolations from here and now to there and then. But, there is also 
the possibility that these challenges to SBBN are providing hints of new physics 
beyond the standard models of particle physics and/or cosmology. Once this 
option is entertained, the possibilities are limited only by the imagination and 
creativity of physicists and cosmologists. Many such models have already been 
proposed and studied. One option is to discard most or all of standard physics 
and start afresh. A more conservative (not in the pejorative sense!) approach 
is to recognize that the standard model (SBBN) does quite well and to look 
for small variations on the same theme. Here, I'll adopt the latter strategy and 
explore one such option which can resolve some, but not all, of the conflicts: a 
non-standard, early universe expansion rate. 

There are many different extensions of the standard model of particle physics 
which result in modifications of the early universe expansion rate (the time - 
temperature relation). For example, additional particles will increase the energy 
density (at fixed temperature), resulting in a faster expansion. In such situations 
it is convenient to relate the extra energy density to that which would have been 
contributed by an additional neutrino with the ordinary weak interactions |19|. 
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Just prior to annihilation, this may be written as 

P' 7Z\N. 
p 43 ■ 



(1) 



Since the expansion rate (the Hubble parameter) depends on the square root of 
the combination of G (Newton's contstant) and the density, the expansion rate 
factor [S) is related to AN^ by, 

= (l + ^)^/^ = (§)^/^. (2) 

Thus, while adding new particles (increasing the energy density) results in a 
speed-up in the expansion rate {AN^ > 0; 5 > 1), changing the effective New- 
ton's constant may either increase or decrease the expansion rate. Another ex- 
ample of new physics which can alter the expansion rate is the late decay of a 
very massive particle which reheats the universe, but to a "low" reheat temper- 
ature 120]. If the reheat temperature is too low (;$ 7 MeV) the neutrinos will 
fail to be fully populated, resulting in Z\Nj^ < and 5 < 1. Finally, in some 
higher dimensional extensions of the standard model of particle physics addi- 
tional terms appear in the 3-f 1 dimensional Friedman equation whose behavior 
mimics that of "radiation" , resulting in an effective ANi, which could be either 
positive or negative PT|. 




Fig. 5. Isoabundance curves for *He (solid) and D (dashed) in the baryon abundance 
(7710) - expansion rate factor (5*) plane. The labels on the '*He curves are for Yp, 
while those on the D curves are for yD = 10^(D/H). The filled circle with error bars 
corresponds to the adopted values of the D and ''He primordial abundances (see the 
text). 
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Thus, a nonstandard expansion rate (S 7^ 1) is a well-motivated, one param- 
eter modification of SBBN which has the potential to resolve some of its chal- 
lenges. A slower expansion would leave more time for neutrons to become protons 
and a lower neutron abundance at BBN would result in a smaller Yp (good!). 
Since "^He is the most sensitive chronometer, the effect on its abundance is most 
significant. However, a modified expansion rate would also affect the predicted 
abundances of the other light nuclides as well. A slowdown will result in more 
destruction of D and '^He and, for 7710 ^ 4, in production of more ^Be (which 
becomes ^Li via electron capture). In Figure are shown the quite accurate ap- 
proximations to the isoabundance curves for D and Yp in the S* — 7710 plane, from 
the recent work of Kneller and myself ^ . Also shown in Fig.Ois the location in 
this plane corresponding to the adopted D and "^He abundances (including their 
uncertainties). Not surprisingly, it is possible to adjust these two parameters (5 
and 77) to fit the relic abundances of these two nuclides. Note, however, that the 
best fit corresponds to a slower than standard expansion rate (S « 0.94 ± 0.03; 
ANi, Ri —0.70 ± 0.35). While this combination of parameters is consistent with 
WMAP (see, e.g., Barger et al. and further references therein), it does not 
resolve the conflict with ^Li. Although a slowdown in the expansion rate does 
have the effect of increasing ^Li, this is compensated by the somewhat lower 
baryon density, which has the opposite effect. The result is that for these choices 
of S and r/, which do resolve the conflicts between D and ^He (and WMAP and 
''He), the predicted primordial abundance of ^Li is still [Li]p « 2.62 ± 0.10. 

4 Summary and Conclusions 

Four light nuclides (D, ■^He, "^He, ^Li) are predicted to emerge from the early 
evolution of the universe in abundances large enough to be observed at present. 
In SBBN there is only one parameter, the baryon density parameter rj, which 
determines the relic abundances of these nuclides. The current observational 
data identifies a range in this parameter of about a factor of 2-3 for which the 
SBBN predictions are in agreement with the primordial abundances inferred 
from current data. This range of t] is also consistent with independent, non- 
BBN estimates 0. Tests of the standard model of cosmology at 20 minutes 
(BBN) and nearly 400 thousand years later (WMAP) agree. While this is a 
great triumph for the standard models of cosmology and of particle physics, the 
agreement is not perfect and, if the uncertainty estimates are taken seriously, 
there are some challenges to this standard model. In this talk, to an audience 
of astronomers, I have emphasized the observational uncertainties in the hope 
of helping to stimulate further observational (and theoretical) work. Will more 
and better data resolve these apparent conflicts? Or, will we be pointed to new 
physics beyond the current standard models? Only time will tell. 
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